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Abstract -The magnetoconductance of graphene nanoribbons with rough zigzag and armchair 
edges is studied by numerical simulations, nanoribbons with sufficiently small bulk disorder show 
a pronounced magnetoconductance minimum at cyclotron radii close to the ribbon width, in 
close analogy to the wire peak observed in conventional semiconductor quantum wires. In zigzag 
nanoribbons, this feature becomes visible only above a threshold amplitude of the edge roughness, 
as a consequence of the reduced current density close to the edges. 



Graphene, i.e. a single layer of graphite, possesses 
an energy dispersion close to the charge neutrality point 
which is equivalent to the Dirac dispersion for masslcss 
fcrmions. [UH] This property has raised tremendous ac- 
tivities in fundamental research in the past few years, 
while the high room temperature electron mobility [3] of 
up to 20, 000 cm 2 /Vs together with excellent mechanical 
properties has raised hopes for novel electronics applica- 
tions. One key element of graphene nanostructures and 
devices are ribbons with submicron width W, also known 
as graphene nanoribbons (GNRs). It is well established 
that in GNRs, edge disorder can contribute significantly 
to the scattering [4HZ] which is governed in wide structures 
by a combination of scattering at charged impurities and 
resonant scattering at short-range defects. [8j[9] Further- 
more, edge roughness has been suggested as the source 
of the transport gap in narrow GNRs around the charge 
neutrality point. [7l fT0TfT2] It would therefore be desirable 
to be able to characterize and quantify the edge roughness 
in GNRs via transport experiments. Conventional quan- 
tum wires made from, e.g., III/V semiconductors, show 
a characteristic magnetoresistance maximum at a perpen- 
dicular magnetic field which corresponds to a cyclotron ra- 
dius r c « Q.hhW , commonly refereed to as wire peak, [13] 
which can be interpreted as an enhanced diffusive scatter- 
ing when the classical electron trajectory grazes at the wire 
edge. |13pi4j Within a quantum picture, the reduced con- 
ductance emerges from a homogenized contribution of the 
occupied wire modes to the overall conductivity. [14j Since 



the strength of this effect depends on the parameters that 
characterize the edge roughness, namely its amplitude and 
correlation length which in turn determine the specularity 
of the edge scattering, [13l[T4] the wire peak is routinely 
used to characterize the edge roughness. [T3lfT5T - H7] We 
are not aware of any observation of a similar structure 
in GNRs, where magnetoresistance peaks in the param- 
eter range of interest are absent. [T5H2Tj This difference 
between GNRs and conventional quantum wires has re- 
mained unexplained. The goal of the present paper is to 
fill this gap and discuss theoretically the edge roughness 
induced magnetoconductance dip (ERID) in GNRs. Our 
results show that this structure should be observable in 
GNRs as well provided the bulk disorder is sufficiently low. 
Furthermore, in zigzag GNRs, the structure is strongly 
suppressed at small edge disorder amplitudes. 



We start with perfect armchair or zigzag nanoribbons 
oriented in x-direction and of width W in y-direction. 
Edge roughness is simulated at both edges by missing 
atoms. It is parameterized by the maximum number 5 
of missing atoms in y-direction from the edges into the 
GNR, and by the correlation length £. The positions of the 
notches in x-dircction arc chosen randomly with a prob- 
ability p, while their individual central depth d c is dis- 
tributed equally between and 8. Then the atoms around 
the defect centers are removed following the Gaussian-type 
shape such that the superimposed notches produce the 
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Fig. 1: (Colour online) Magnetoconductance of armchair GNRs 
(width 50 nm, length 220 nm) as a function of the edge rough- 
ness parameters S (a) and £ (b). Each trace is averaged over 
500 edge configurations. The dashed lines show the number 
of occupied subbands. Inset in (a): typical edge profiles for 
specific parameter values. Inset in (b): amplitude of the ERID 
as a function of £. 



number of missing atoms at edge site Xj 

d(x,) = ^e X p(-^^) (1) 

Here, N c — p ■ N at0 m is the number of notch centers with 
Natom being the number of carbon atoms along the edges. 
The roughness patterns on both edges are uncorrelated. 
Typical structures are shown for one edge in the inset of 
fig.Hl(a). 

Besides the edge roughness, we also study the ERID ef- 
fect under the influence of bulk disorder which is modeled 



by the widely used random Gaussian potential [251150] 

V^Y^U.exJ-^fL) (2) 

i' = l ^ ' 

The potential heights Ui are uniformly distributed in the 
range Ui G [—A, A] where A denotes the maximum poten- 
tial height. The potential decay depends on the site-site 
distance |r.; — | and the correlation length A. We refer to 
the case A = a (appropriate for the absorbed neutral im- 
purities) as short-range scattering, and to the case A>s 
(appropriate for the remote charged impurities) as long- 
range scattering, with a = 0.142nm being the C-C bond 
length. 

Our calculations are based on the tight-binding Hamil- 
tonian of GNRs 

with the hopping energy = t = — 2.7eV between 
the nearest neighboring carbon atoms at zero magnetic 
field. The perpendicular magnetic field is incorporated via 

ie/h f J A. d\ 

Peierls' substitution — > i^e J > with the vector 
potential A = {—By, 0, 0) in the Landau gauge. 

The conductance of GNRs is calculated with Landauer- 
Buttikcr formula G = 2e 2 /hT(Ep), where the transmis- 
sion T(Ep) is computed with the aid of the recursive 
Green's function technique [B]. Throughout this work, the 
Fermi energy Ep = O.lt corresponding to an electron den- 
sity of n = 5.3 x 10 16 m~ 2 is used. The Fermi wavelength 
Xf is obtained from the linear band relation Ep = TivfUf 
and Xp = 2ir/kp with the Fermi velocity vp = 3ta/2h. 
The cyclotron radius is given by r c = Tikp/eB. 

Fig. [U shows the results for the magnetoconductance 
G(B) of armchair GNRs of width W = 50 nm and length 
L = 220 nm for various values of the edge roughness pa- 
rameters and in the absence of bulk disorder. As the 
magnetic field is increased from zero, G(B) decreases and 
forms a broad minimum, the ERID, at B m i n « 5 T corre- 
sponding to W/r c ~ 0.8. By varying W in the simulations, 
we find that B min is nearly inversely proportional to the 
GNR width W (not shown). This structure is the ana- 
logue of the so-called wire peak observed in conventional 
semiconductor quantum wires. |13j The amplitude of the 
ERID depends only weakly on the edge roughness ampli- 
tude S, Fig. [T{a), and decreases as the correlation length 
of the roughness £ is increased, Fig. [TJb). Outside the 
dip, a substantial increase of the conductance follows as 
B is increased, corresponding to the positive magnetocon- 
ductance observed in recent experiments. [T5JI2D] As the 
magnetic field is increased further beyond « 12 T, oscilla- 
tory structures appear which reflect the magnctodepletion 
of the GNR modes and are of no further interest here. We 
note that the effective widths of GNRs are reduced by 
rough edges such that their signature in the magnetocon- 
ductance is slightly displaced with respect to the mag- 
netodeplction of modes in defect-free GNRs. The ERID 
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Fig. 2: (Colour online) Magnetoconductance of a GNR with 
zigzag edges and otherwise identical parameters to those in Fig. 
1, for various edge disorder amplitudes and a fixed correlation 
length of £ = 2a (a), where the dashed line denotes the number 
of occupied modes, and for different correlation lengths at a 
disorder amplitude of S = 10 (b). The inset in (b) shows the 
cross-sectional current density for zigzag and armchair GNRs 
with perfect edges. 



occurs in an interval where the number of occupied modes 
does not change much, which indicates that the structure 
is not related to the diamagnetic shift of the modes. This 
behavior is in qualitative agreement with the results of 
the quantum treatment of the wire peak in conventional 
quantum wires by Akera and Ando. [14] 

Fig. [2]shows the corresponding results for zigzag GNRs. 
Again, the ERID structures are observed at cyclotron radii 
somewhat larger than W, but distinct differences to arm- 
chair GNRs can be identified. First of all, at B = 0, 
the conductance suppression by the edge disorder is less 
pronounced as compared to armchair GNRs. This repre- 
sents the well-known reduced sensitivity to edge disorder 
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of zigzag GNRs compared to armchair GNRs. [3]. In ad- 
dition, at small disorder amplitudes, e.g. for (5 = 3, up- 
permost trace in Fig. [Ha), the ERID is absent and the 
magnetoconductance reflects just the depopulation of the 
modes. As S increases to 6, the ERID becomes visible, 
and its amplitude approaches that one observed in arm- 
chair GNRs as 5 is increased. As £ is increased at an 
amplitude of 6 = 10, Fig. [2] (b), the ERID amplitude 
develops nonmonotonically and shows a minimum around 
£ = 4a. For larger correlation lengths, the ERID develops 
an asymmetric shape. 

Apparently, compared to armchair GNRs, the magne- 
toconductance in zigzag GNRs is less sensitive to edge 
disorder for small disorder amplitudes. In order to inter- 
pret this difference, we take a look at the current density 
distribution j x (y) in GNRs with perfect edges, see the in- 
set in Fig. [5] (b). The current density of armchair GNRs 
oscillates as a function of the site index (red trace), which 
is a consequence of the alternating character of the wave 
functions. Averaging over 4 adjacent atoms (blue trace) 
however shows that the current density is almost homo- 
geneously distributed over the cross section. In zigzag 
GNRs, however, the current density is strongly suppressed 
in an interval extending about « 5 ran from the edge into 
the bulk. This current density suppression, which has been 
reported earlier, [35] suggests that zigzag GNRs are less 
sensitive to edge disorder, in particular for disorder ampli- 
tudes below s=s 5 nm for the Fermi energy considered here, 
which is in good qualitative agreement with the magneto- 
conductance simulations. Furthermore, it has been shown 
that in the single-mode regime in zero magnetic field, the 
conductance of the zigzag GNRs remains practically unaf- 
fected by the edge disorder [3] . In this case, the propagat- 
ing channels in the zigzag GNRs are edge states flowing 
on the opposite boundaries [281129] such that the backscat- 
tcring is suppressed. We believe that this mechanism is 
not relevant in our study since we are in the multimode 
regime where the special character of the edge state should 
be insignificant. It would thus be interesting to study the 
magnetoconductance very in close proximity to the charge 
neutrality point. Our numerical method is an impractical 
tool for such a study due to the onset of Coulomb block- 
ade effects as the Fermi energy is reduced. [57] We did, 
however, vary the number of occupied subbands at fixed 
Fermi energy by changing W and found that the ERID 
remains located close to r c = W. 

To discuss implications regarding the experimental ob- 
servability of the ERID, we first study the effect of bulk 
disorder. The edge disorder of an armchair GNR is fixed 
to the parameters 5 — 10 and £ = 2a as used in Fig. 
HI Fig- El (a) shows the magnetoconductance for various 
values of the bulk disorder potential amplitudes A, for 
short-range disorder (A = a) and with the impurity den- 
sity riimp = 8x 10 12 cm~ 2 . As A gets larger, the absolute 
amplitude of the ERID gets suppressed while its position 
remains constant. In addition, a small conductance dip 
emerges at B = as the bulk disorder is increased. We in- 
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Fig. 3: (Colour on-line) (a) Effect of bulk disorder on the 
edge roughness induced magnetoconductance dip in armchair 
GNRs, averaged over 500 configurations. The width and length 
of GNR are 50 nm and 220 nm, respectively. The edge disorder 
parameters are 8 = 10 and £ = 2a, while the bulk disorder is 
varied by tuning impurity center height A. The inset shows the 
ERID amplitude as a function of A. (b) Length dependence of 
the magnetoconductance. 

tcrprct this as the onset of bulk-disorder induced weak lo- 
calization. This behavior remains the same for long range 
potentials (A = 32a, not shown). 

Second, the length dependence of the ERID is studied. 
As the length L of the GNR is increased from 100 nm 
to 200 nm, the ERID amplitude increases and tends to 
saturate as L is increased further. We take this as an 
indication that the system goes from the ballistic into the 
diffusive regime around a length of 200 nm. We note that 
in all cases studied here, the localization length ti oc is 
larger than L. For the case shown in Fig. [3](b), li oc = 2fim 
is found by fitting the scaling law ln(l + l/g) = L/£i oc with 
g denoting the dimensionless conductance g = G/(e 2 /irh). 



[HUGH] For GNRs with L > l\ oc , the conductance would 
be exponentially suppressed. 

These investigations allow us to draw some conclusions 
regarding the observability of the ERID in GNRs in exper- 
iments and concerning the suitability of this effect to char- 
acterize the edge roughness. It has emerged that the ERID 
can be observed in GNRs only if the bulk disorder level is 
sufficiently low. This may not be the case in present ex- 
perimental implementations. Typical experimental edge 
roughness amplitudes arc in the range of 1-10 missing 
atoms, which is the range where the wire peak should 
be observable in armchair GNRs provided the correlation 
length of the edge roughness remains small, i.e., of the or- 
der of the Fermi wavelength or below. In zigzag GNRs, on 
the other hand, the ERID will be only visible if the edge 
roughness amplitude is not small compared to the interval 
of suppressed current density at the GNR edges, i.e. S > 6 
for the parameter range studied here. Like in conventional 
quantum wires, scattering at the edges gets less diffusive 
if the roughness becomes smoother which suppresses the 
ERID. Furthermore, the GNR lengths should be smaller 
than the localization length but well within the diffusive 
regime for maximum visibility. We note that there is no 
dephasing mechanism in our model. In real samples with 
finite dephasing time, the conductance around B=0 will 
increase due to the suppression of weak localization. 

In summary, we have studied the magnetotransport 
properties of both armchair and zigzag GNRs in the pres- 
ence of edge disorder. In both types of GNRs, an edge 
roughness - induced conductance minimum is found at 
cyclotron radii corresponding roughly to the GNR width, 
provided the edge disorder has a correlation length not 
much larger than the Fermi wavelength and the bulk dis- 
order is sufficiently low. The effect is most visible in the 
diffusive regime. This magnetoconductance structure is 
similar in character to that one observed at conventional, 
rough quantum wires. Furthermore, in zigzag GNRs, it 
gets suppressed as the edge disorder amplitude drops be- 
low the spatial extension of the reduced current density. 
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